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Abstract 



Graphene is considered one of the most promising materials for future electronic. However, in its pristine 
form graphene is a gapless material, which imposes limitations to its use in some electronic applications. 
In order to solve this problem many approaches have been tried, such as, physical and chemical functional- 
izations. These processes compromise some of the desirable graphene properties. In this work, based on ab 
initio quantum molecular dynamics, we showed that a two-dimensional carbon allotrope, named bipheny- 
lene carbon (BPC) can be obtained from selective dehydrogenation of porous graphene. BPC presents a 
nonzero bandgap and well-delocalized frontier orbitals. Synthetic routes to BPC are also addressed. 



1 Introduction 

In the last decades the successive discoveries of new 
carbon-based materials have opened a new era in 
materials science. Examples of these discoveries are 
fullerenes [16 , carbon nanotubes [14] . and more re- 
cently, graphene [22] |2T], 

Graphene is a two-dimensional array of hexago- 
nal units of sp 2 bonded carbon atoms (Fig. la). 
Graphene presents very unusual and interesting elec- 
tronic and mechanical properties [22j [21] . Because of 
these special properties graphene is considered one of 
the most promising materials for future electronics. 
However, in its pristine form graphene is a gapless 
semiconductor, as shown in Fig. 2a. This poses seri- 
ous limitations to its use in some electronic applica- 
tions, such as, some kind of transistors [3T] . 

Many approaches have been tried in order to create 
a gap in graphene-like materials. The most common 
strategies use chemical and physical methods, such 
as, oxidation [25], hydrogenation [28, 29j[10][T2] and 
fluorination [6j[20l[T7]. However, the controlled syn- 
thesis of large structures and/or at large scale have 
been proved to be difficult. More important, the 
desirable electronic graphene properties are partially 
compromised in such approaches. 



Another approach has been trying to obtain intrin- 
sically hydrogenated structures, such as, the so-called 
porous graphene (PG) (Fig. lb) 01 El El H3 HB], 
whose synthesis has been recently achieved [4]. But 
again, the obtained structures present some of the 
same problems of chemically /physically functional- 
ized graphene, such as, excessive large bandgap value 
and flat (low mobility) electronic bands (Fig. 2b). 

An ideal structure would be an allotrope carbon 
form with an intrinsic good bandgap value and elec- 
tronic bands with good dispersion (electron dereal- 
ization and high charge mobility). In theory, struc- 
tures satisfying these conditions do exist, as the so- 
called biphenylene carbon (BPC) (Figs, lc and 2c) 
[2] [n]. Molecular fragments in linear, zigzag and 
several other forms, have already been synthesized 
[26] . They are considered as potential precursors for 
fullerenes, bowls, cyclophanes, etc. [2j[TTJ[26]. How- 
ever, the synthesis of large BPC fragments remains 
elusive. 

In this work, based on ab initio molecular dynam- 
ics simulations, we show that selective dehydrogena- 
tion of porous graphene leads to the spontaneous in- 
terconversion to BPC structures. Possible synthetic 
route approaches to achieve this interconversion are 
also discussed. That this interconversion was possi- 
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Figure 1: (Color online) Structural models considered 
in the present work, (a) Graphene, with two atoms 
in the unit cell, (b) porous graphene with 18 atoms in 
the unit cell and (c) BPC with 12 atoms in the unit 
cell. Unit cell is highlighted in each case. 

ble was found out by exploratory investigations, as 
discussed below. 



2 Methodology 

The geometric and electronic aspects of the struc- 
tures shown in Fig. 1, as well as, the interconversion 
processes from PG to BPC, were investigated in the 
framework of density functional theory (DFT). Ex- 
change and correlation terms were considered within 
the generalized gradient approximation (GGA) with 
a BLYP functional [3 , with a double numerical plus 
polarization basis set, as implemented in DMol 3 code 
[7J [8]. In all the calculations a non-relativistic and 
all-electron treatment was used. The parameter cri- 
teria for the tolerances of energy, force, displacement 
and SCF convergence criteria were 2.72 x 10 -4 eV, 




Figure 2: (color online) Band structure (in eV) for 
the (a) graphene, (b) porous graphene and (c) BPC. 
Dashed red lines indicate the Fermi level. Please no- 
tice that DFTB+ convention to locate the Fermi level 
is at the value of the highest occupied state. 

5.44 x 1(T 2 eV • A~\ 5.0 x 1(T 3 A and 1.0 x 10~ 6 , 
respectively. 

It is well-known that although DFT methods can 
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Table 1: DMol 3 geometrical data. Labels according 
to Fig. ^ 

GRAPHENE 



Distances (A) 


Angles (°) 


Ci-C 2 1.39 


C1-C2-C3 119.96 



POROUS GRAPHENE 



Distances (A) 


Angles (°) 


C1-C2 1.41 


Ci- 


-C 5 -C 6 121.56 


C3-C4 1.41 


c 2 - 


-C1-C5 122.68 


C 5 -C 6 1.50 


Cr 


-C2-C3 117.77 






BPC 


Distances (A) 




Angles (°) 


C1-C2 1.476 


Ci 


-C2-C3 89.826 


C1-C4 1.485 


c 2 


-C3-C4 90.174 


C1-C5 1.365 


c 2 


-C1-C5 119.846 



reliably describe the geometrical features, in general 
some electronic properties, as the bandgap values, are 
underestimated [27]. For these reasons, for the anal- 
ysis of the electronic band structure calculations we 
used a DFT based tight-binding method (DFTB+) 
[24j [1], which has been proved to reliably describe 
the electronic properties of carbon-based materials. 
The DFTB+ calculations were carried out using the 
optimized geometries from the DMol 3 calculations. 
The DFTB+ non-diagonal matrix (overlap matrix 
and Hamilton matrix) elements are calculated in a 
two-center approximation. They are distant depen- 
dent and considered up to about 10 atomic units. 

Initially, we carried out DMol 3 fully (unit cell pa- 
rameters and atomic positions were allowed to vary) 
geometric optimizations. The most relevant geomet- 
rical data are displayed in Tables I and II. See also 
supplementary materials. 

3 Results 

Our results show that the PG optimized geometry is 
characterized by a (7222(^1) symmetry group. The 
hexagonal geometry is preserved, with bond-length 
in the rings and intra-rings of 1.40 and 1.50 A, re- 
spectively. The obtained geometry is in good agree- 
ment with the experimental data reported [4], theory 
and experiment estimate the lattice parameter to be 
about 7.4 A. The obtained bandgap value of 3.3 eV 
is also in good agreement with previous theoretical 



Table 2: Lattice parameters for graphene, porous 

graphene and biphenylene carbon (BPC), respec- 
tively. 

a (A) b(A) 

Graphene 2.40 2.40 

Porous Graphene 7.52 7.53 

BPC 6.78 6.69 



calculations (3.2 eV) [9]. 

For the BPC layer, a stable conformation was ob- 
tained and satisfying the topological conditions of the 
theoretically proposed structures [2j [11] . In relation 
to PG, there is a significant lattice parameter con- 
traction, from 7.5 to 6.8 A, respectively (Table II). 
The bond-lengths in the hexagons preserve the pat- 
tern of alternating double and single bonds (1.48 and 
1.36 A, respectively), while the square structure (cy- 
clobutadiene) is almost a perfect square (1.48 A) (see 
Fig. 1 and Table I). 

In Fig. 2 we present the band structure results 
obtained from DFTB+ calculations. As expected, 
graphene presents a zero bandgap value, while PG 
and BPC have values of 3.3 and 0.8 eV, respectively. 
However bandgap values and good dispersion of the 
frontier bands are not warranty of good conductors. 
Besides these aspects another important character- 
istic which plays an important role in defining the 
electronic (conductivity) mobility of the material is 
the degree of electronic derealization of the frontier 
crystalline orbit als. The electronic analysis we car- 
ried out here for ideal BPC structures are for neutral 
forms (no free carriers). In order to create these car- 
riers is necessary, as usual, to dope the material. In 
ref [16] the doping of porous graphene has been ad- 
dressed. The same principles can be used to dope 
BPC structures. 

In Fig. 3 we present the frontier orbitals HOCO 
(highest occupied crystalline orbital) and LUCO 
(lowest unoccupied crystalline orbital) for PG and 
BPC. As we can see from Fig. 3 BPC presents more 
delocalized orbitals than PG. Particularly interesting 
is the BPC HOCO, which is well delocalized over the 
whole network, suggestive of a good conductor struc- 
ture. Thus, in principle, BPC presents itself as an 
ideal structure for many electronic applications; good 
intrinsic bandgap value, bands with good dispersion 
and delocalized frontier orbitals. Also, we calculated 
the values of the effective masses (m*/m e =0.26 and 
m* jvfih = 0.33, for the conduction and valence ef- 
fective masses, respectively). These values are con- 
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Figure 3: (Color online)DFTB+ frontier orbitals for 
hydrogenated and dehydrogenated porous graphene. 
(a) and (b) HOCO and LUCO for porous graphene. 
(c) and (d) HOCO and LUCO for dehydrogenated 
porous graphene. 



sistent with the expected to a good conductor. How- 
ever, as mentioned before, the synthesis of large BPC 
fragments from using chemical methods remains elu- 
sive. Through exploratory investigations we found 
out that selective hydrogen removal from PG leads 
to a spontaneous interconversion to BPC. 

We have carried out a b initio quantum molecular 
dynamics (DMol 3 ), NVT ensemble, at different tem- 
peratures: (just geometric optimizations), 300 and 
600 K. We started from the optimized PG geometry, 
then we removed the hydrogen atoms and let the sys- 
tem to freely evolve in time (lattice parameter values 
and atomic positions free to vary). For all the inves- 
tigated cases we observed a spontaneous interconver- 
sion from dehydrogenated PG to BPC. The obtained 
BPC structures are thermally stable (at least up to 
600 K). 

In Fig. 4 we present results for the calculations car- 
ried out at 300 K. In Figs. 4a-4d we show snapshots 
from the molecular dynamics simulations at succes- 
sive time steps. We can see that the dehydrogenated 
PG undergoes structural rearrangements, mainly ring 
rotations (about 30 degrees) coupled to a lattice pa- 
rameter reduction leading to the formation of a cy- 
clobutadiene motif (highlighted atoms in Figs. 4a- 
4d), and consequently to the BPC formation. The 
whole process can be better visualized in the videoOl 




Figure 4: DMol 3 molecular dynamics snapshots and 
total energy values as a function of time, (a) Com- 
pletely dehydrogenated porous graphene on the ini- 
tial stage. The highlighted carbons in red lead to the 
formation of a cyclobutadiene. (b) and (c) Inter- 
mediate stages. Notice that adjacent rings rotate in 
opposite directions, (d) Final stage after a complete 
rotation of 30° in each ring, (e) Time evolution for 
the total energy during the QMD simulation. Letter 
labels represent the energy of each snapshot. 



of the supplementary materials. 

In Fig. 4e we present the total energy values as a 
function of time of simulation. As we can see, start- 
ing from dehydrogenated PG the system continuously 
evolve to more stable configurations reaching a well- 
defined minimum, which is associated with the BPC 
formation. These results can be explained in terms 
of the relative stability of graphene, PG and PBC. 
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The total energy values per carbon atoms (in re- 
lation to graphene) of dehydrogenated PG and BPC 
are 1.26 and 0.63 eV, respectively. In this sense the 
removal of hydrogens from PG inverts the stability 
order in relation to BPC, BPC being now 0.63 eV 
per carbon atom more stable. Due to the similar 
topology the interconversion easily occurs, since it 
requires only ring rotations and the creation of new 
bonds forming the cyclobutadiene motif. These re- 
sults strongly support that BPC can be obtained from 
PG by just selective dehydrogenation. 

PG selective dehydrogenation is within our 
present-day synthesis capabilities. Recently, signifi- 
cant experimental advances have been produced in 
selective dehydrogenation of hydrocarbons [23j [30j 
\19\ [T3] . Fullerenes [23], nanographene flakes [30] , 
and even two-dimensional networks [I3j [19] similar to 
BPC have been already achieved. See more details in 
the supplementary materials. Obtaining BPC from 
DPG using these techniques is perfectly feasible. We 
hope the present results will stimulate further works 
along these lines. 

4 Acknowledgements 

This work was supported in part by the Brazilian 
Agencies CNPq, CAPES and FAPESP. 

5 Supporting Information 

Brief discussion about simulations and MD videos. 
This information is available free of charge via Inter- 
net at http://pubs.acs.org 

References 

[1] B. Aradi, B. Hourahine, and T. Frauenheim. J. 
Phys. Chem. A, 111:5678, 2007. 

[2] R. Baughman, H. Eckhardt, and M. Kertesz. J. 
Chem. Phys., 87:6687, 1987. 

[3] A. D. Becke. J. Chem. Phys., 88:2547, 1988. 

[4] M. Bieri and et al. Chem. Commun., 6919, 2009. 

[5] S. Blankerbur, M. Bieri, R. Fasel, K. Mullen, 
C. A. Pignedoli, and D. Passerone. Small, 
6:2266, 2010. 

[6] S. H. Cheng and et al. Phys. Rev. B, 81:205435, 
2010. 



[7] B. Delley. J. Chem. Phys., 92:508, 1990. 

[8] B. Delley. J. Chem. Phys., 113:7756, 2000. 

[9] A. J. Du, Z. H. Zhu, and S. C. Smith. J. Amer. 
Chem. Soc., 132:2876, 2010. 

[10] D. Elias and et al. Science, 323:610, 2009. 

[11] A. N. Enyashin and A. L. Ivanovskii. Phys. St. 
Sol. B, 1, 2011. 

[12] M.Z. S. Flores, P. A. S. Autreto, S. B. Legoas, 
and D. S. Galvao. Nanotechnology, 20:465704, 
2009. 

[13] R. Gutzler and et al. Chem. Commun, page 
4456, 2009. 

[14] S. Iijima. Nature, 354:56-58, 1991. 

[15] D. E. Jiang, V. R. Cooper, and Z. F. Chen. Nano 
Lett, 9:4019, 2009. 

[16] H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. 
Curl, and R. E. Smalley. Ceo- Buckminster- 
fullerene. Nature, 318:162-163, 1985. 

[17] O. Leenaerts, H. Peelaers, AD Hernandez- 
Nieves, B. Partoens, and FM Peeters. First- 
principles investigation of graphene fluoride and 
graphane. Physical Review B, 82(19):195436, 
2010. 

[18] Y. F. Li, Z. Zhou, P. W. Shen, and Z. F. Chen. 
Chem. Commun, 46:3672, 2010. 

[19] H. Liang and et al. Coord. Chem. Rev., 253:2959, 
2009. 

[20] R. R. Nair and et al. Small, 6:2877, 2010. 

[21] A. H. Castro Neto, F. Guinea, N. M. R. Peres, 
K. S. Novoselov, and A. K. Geim. Rev. Mod. 
Phys., 81:109, 2009. 

[22] K. Novoselov and et al. Science, 306:666, 2004. 

[23] G. Otero and et al. Nature, 454:865, 2009. 

[24] D. Porezag, T. Frauenheim, T. K. Kohler, 
G. Seifert, and R. Kaschner. Phys. Rev. B, 
51:12947, 1995. 

[25] R. Ruoff. Nature Nanotech, 3:10, 2008. 

[26] J. M. Schulman and R. L. Disch. J. Phys. Chem. 
A, 111:10010, 2007. 



5 



[27] David Sholl and Janice A. Steckel, editors. Den- 
sity Functional Theory: A Practical Introduc- 
tion. Willey-Inter science, New York, 2009. 

[28] M.H.F. Sluiter and Y. Kawazoe. Cluster expan- 
sion method for adsorption: Application to hy- 
drogen chemisorption on graphene. Physical Re- 
view B, 68(8):085410, 2003. 

[29] J. O. Sofo, A. S. Chaudhari, and G. D. Barber. 
Phys. Rev. B, 75:153401, 2007. 

[30] M. Treier and et al. Nature Chem., 3:61, 2011. 

[31] F. Withers, M. Dubois, and A. K. Savchenko. 
Phys. Rev. B, 82:73403, 2010. 



6 



